Quantum Anomalous Hall Effect in Flat Band Ferromagnet 



An Zhao and Shun-Qing Shen^ 

'Department of Physics and Center of Computational and Theoretical Physics, 
The University of Hong Kong, Pokfulam Road, Hong Kong 

We proposed a theory of quantum anomalous HaU effect in a flat-band ferromagnet on a two- 
dimensional (2D) decorated lattice with spin-orbit coupling. Free electrons on the lattice have 
dispersionless flat bands, and the ground state is highly degenerate when each lattice site is occu- 
pied averagely by one electron, i.e., the system is at half filling. The on-site Coulomb interaction 
can remove the degeneracy and give rise to the ferrimagnetism, which is the coexistence of the 
ferromagnetic and antiferromagnetic long-range orders. On the other hand the spin-orbit coupling 
makes the band structure topologically non-trivial, and produces the quantum spin Hall effect with 
a pair of helical edge states around the system boundary. Based on the rigorous results for the 
Hubbard model, we found that the Coulomb interaction can provide an effective staggered potential 
and turn the quantum spin Hall phase into a quantum anomalous Hall phase. 

PACS numbers: 71.10.Fd, 71.70.Ej, 75.50.Ee 



I. INTRODUCTION 

Quantum anomalous Hall effect (QAHE) is a quantum 
mechanical version of the Hall effect in a ferromagnet in 
absence of an external magnetic field or Landau levels. 
Different from the quantum Hall effect in a strong mag- 
netic field, it originates from the topological properties 
of band structure in solid. Usually the anomalous Hall 
effect occurs in ferromagnetic metals due to the intrin- 
sic spin-orbit coupling or extrinsic spin-orbit scattering ii 
It has been realized that the Hall conductance can be 
expressed as the summation of the Berry curvature in 
momentum space over all occupied states of electrons. 
This makes it possible to realize the quantum Hall effect 
even without the Landau levels.^ A simple picture for 
QAHE was proposed in a two-dimensional electron gas 
with strong spin-orbit coupling.— Electron localization in 
disordered systems also provides an alternative approach 
to realize this effectj^i^ The discovery of the quantum 
spin Hall effect (QSHE) and topological insulators^"— 
stimulated extensive interest to search QAHE in realistic 
systems. As a result of time-reversal symmetry break- 
ing in topological insulators, this effect was predicted 
in a magnetically doped Hgi_j^MnyTe quantum welli^ 
and Cr or Fe doped topological insulators in which 
the presence of magnetization suppresses one of helical 
edge states in QSHE, and preserves a chiral edge state 
for the quantum Hall effect. Very recently HgCr2Se4 is 
proposed to be a Chern semimetal and exhibits QAHE 
in quantum- well structure>i^ Interaction-driven topologi- 
cally non-trivial Mott insulating phase displaying QAHE 
or QSHE has also attracted much research interest. ^^"^^ 

While in search of ferromagnetism in diluted magnetic 
semiconductors, rigorous models of ferromagnetis m^^i^° 
or ferrimagnetis m^ in strongly correlated electron sys- 
tems provides an alternative context to realize ferromag- 
netism. A common feature of these examples is the flat or 
almost flat band of electrons. The Coulomb interaction 
may remove high degeneracy of electrons in the band, 
leading to the ferromagnetism according to the Stoner 



criteria. Recently it was predicted that the ferromag- 
netism may occur in the Hubbard model with topological 
non-trivial bands. 

In this paper we propose a theory of QAHE in a flat 
band ferromagnet with spin-orbit coupling. We start 
with a two-dimensional decorated lattice model, which 
has a pair of flat bands in the middle of the energy spec- 
trum. The inclusion of spin-orbit coupling makes the 
energy bands topologically non-trivial, and gives rise to 
QSHE. Based on the rigorous results for the Hubbard 
model, the Coulomb interaction may induce the ferri- 
magnetism in the ground state, when the middle flat 
bands are half filled. It provides an ideal way to realize a 
magnetic staggered field. The staggered field can modu- 
late the topological numbers of electron bands by closing 
and reopening the energy gapi^ Different configurations 
of topological invariants assigned to the electron bands 
can then be obtained and will give different topological 
phases. Based on a self-consistent mean field calculation, 
we present the phase diagram of the ground state. We 
find that QAHE with nonzero Chern number can be re- 
alized in a ferromagnet due to the Coulomb interaction. 

II. A DECORATED LATTICE MODEL 

We begin with the tight-binding Hamiltonian on a two- 
dimensional decorated lattice [see Fig. [Ha)], 

n = no + nso, (i) 

where the spin-independent hopping term is given by 
Ho ^t^cjcj, 

Ci = (ci,t: Ci.j,)"^ and c|^^^ are the annihilation and cre- 
ation operators of electron with spin t (i) on site i. {i, j) 
means the summation over the nearest-neighbor sites, t 
is the hopping amplitude. The spin-orbit coupling term 
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FIG. 1: (Color online) (a) A two-dimensional decorated lat- 
tice. Dashed lines represent the hopping with spin-orbit cou- 
pling. The light (green) and dark (black) dots denote the sub- 
lattices A and B, respectively, (b) The energy dispersion for 
H(k), which has three pairs of doubly degenerate bands, and 
the gaps open in the presence of spin-orbit coupling (A 7^ 0). 



a. Since H-^ and are time-reversal partners, we 
have C(n,t) = — C(n,4,) for each time-reversal pairs of 
bands, which are degenerate due to time-reversal sym- 
metry. Therefore, the total Chern number is zero. In the 
presence of spin-orbit coupling, the Chern numbers of the 
three bands in {Hi) are {77,0,-77} ({—7/, 0,7/}) from 
top to bottom, with 77 — sgn(A). The nonzero difference 
between C(77, f) and C{n, ].) is equivalent to a non-trivial 
Z2 index, which can also be calculated explicitlyjii^Si"— 
When the Fermi level is located in the gap, the non-trivial 
Z2 index for the filled pairs of bands indicates QSHEi^ 



III. EFFECT OF THE STAGGERED 
POTENTIAL 



has the form 



c] [{dij xdii) ■ cr] ci, 



which gives a spin-dependent hopping between the next- 
nearest-neighbor sites with hopping amplitude A [shown 
by the dash lines in Fig. [IJa)]. Here the lattice is divided 
into two sublattices, A and B, shown by the light and 
dark dots in Fig. (Ha), respectively, j and I denote the 
adjacent sites of site i, and dij is the unit vector along 
the direction from site i to site j, and cr are the Pauli 
matrices. 

We choose the sites 1, 2, and 3 in Fig. [Ua) as the 
unit cell. Since the z component of spin ct^ commutes 
with this Hamiltonian, the Hamiltonian has a block- 
diagonalized form in momentum space after the Fourier 
transformation. 



kcr 



where vj/^^ ^ (ci,ko-, C2,kcr, ca^ko-)'^, cr =t, i denotes dif- 
ferent spins, and 



^t(k) 




2t cos kx 




2t cos ky 4iA sin kx sin 



2t cos ky 
— 4iA sin kx sin 




ky 



Hi(k) = iJ|(— k) is the time-reversal partner of The 
Brillouin zone spans over < fc^; < tt and < ky < n. 
This Hamiltonian preserves time-reversal symmetry, i.e., 
eiJ(k)e-i = H{-k), where 6 = iayK and K is the 
complex conjugate operator. In this spin-1/2 system, the 
time-reversal operator satisfies 8^ = —1. When A 7^ 0, 
in each Ha^, all the three bands are well separated and 
can be characterized by the Chern number— 



<Cin,a) = :^ f d'k [Vx A",(k)] 
Jbz 



(3) 



where ^^^(k) = -i«(k)|Vk|<(k)) and <(k) is the 
Bloch function for the 77th band of electrons with spin 



For an intuitive illustration on how QAHE arises in 
this system, we first introduce a spin-dependent stag- 
gered potential term 



and a spin-independent staggered potential term 



Uc = -V, 



E 



clc,. 



(4) 



(5) 



respectively, where the summations run over the sublat- 
tice sites. In momentum space, the Hamiltonian -ff(k) 
has a block-diagonalized form. 



Hik) 



H[(k) 







(6) 



(2) with H\ = + S+A and H[ = H^- S^A, where A = 



diag(— 1,1,1) and S± — Vs ± Vc- When Vs ^ 0, the 
time-reversal symmetry is broken and the degeneracy of 
the time-reversal pair of bands is removed. In the block- 
diagonalized form, we may say that H-^ feels a staggered 
potential of amplitude S+ while feels —S-. These two 
parts of the Hamiltonian can be investigated separately. 

We notice that a staggered potential may change the 
Chern numbers of the bands of 77| and by closing 
and reopening the band gap in a band inversion. For 
H!^, the bands cross only at the point ko = (7r/2,7r/2) 
when t ^ 0, X ^ 0. The eigenvalues at this point are 
{— 5*+, — 4A + 5'+,4A + S+}, respectively. As a result, 
with increasing 5"+ from zero, a band crossing happens 
at S+ = +2X or S+ = — 2A. For example near S+ — 2A, 
we may obtain an effective Hamiltonian near the point 
ko 

HcB = im + \q^)a^ - \/2t{qxax - qyOy) + X{q^ - 2), 



where q = k — ko and 



^4 



2A. From this two- 



band massive Dirac model it is known that the topologi- 
cal quantum phase transition occurs when the sign of m 
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S =0 — ► S =21. — ► S >2\ 



FIG. 2: (Color online) The band dispersions for and 
along G [0, tt], at ky — n/2 . The horizontal axis is kx and 
the vertical axis is E/t. We take the parameter A > 0. and 
±1 indicate the Chern number of each band. From the left to 
right, < S± < 2A, S± = 2A, and S± > 2A. H_i with 5'- < 
has the identical bands structure as Hf with S+ = —S-, but 
opposite Chern numbers. 

or A changes. The Chern number of the lower band is 
given by^^ (the upper one has a sign change) 

Cl [sgn(m) + sgn(A)] . 

Thus the sign change of m indicates that the Chern num- 
ber changes from to 1 or 1 to 0. Fig. [2] depicts the 
band structure for and with the Chern numbers 
also denoted for each band. There are three different 
cases: (a) when < 5*+ < 2A, the Chern numbers of the 
bands in iJ| are {1,0,-1} from top to bottom; (b) when 
S+ — 2A, the two lower bands touch at kg — (7r/2, 7r/2) 
and the Chern numbers disappear as the two bands are 
not well separated; (c) when S+ > 2A the band gap re- 
opens and the Chern numbers become {1,-1,0} after 
the two inverted bands exchange their Chern numbers. 
Similarly, for the Chern numbers of the three bands 
arc {-1,0,1} when < 5_ < 2A and {0,-1,1} when 
S- > 2A. It is noted that the band structure and Chern 
numbers of fl| with S- and A are identical to those of 
H!^ with 5+ = — 5_ and —A. 

According to the Chern numbers, the edge-bulk corre- 
spondence tells that the edge states in a sample of strip 
geometry depend on parameters. In Fig. [3J we only 
present the edge state spectra of When < 5*+ < 2A, 
the edge state spectra connect the upper and lower bands. 
Since the Chern number of the middle band is zero, this 
topologically trivial band only distorts the edge state 
spectra, but does not affect the existence of the edge 
states. The total Chern number is still 1 and the system 
is topologically non-trivial when the middle and lower 
bands are fully filled. When 5+ > 2A, the middle band 
becomes topologically non-trivial, and the lower band be- 



comes trivial. The edge state spectra only connect the 
middle and the upper bands. 

Now we can present the evolution of edge states in the 
total Hamiltonian H!^ + _f/|. Five cases are listed in Fig. 
m Without loss of generality, we take Vs > and Vc> 0- 
In this case we have 5*+ > \S^\. When time- reversal sym- 
metry is broken, an energy gap can open between a time- 
reversal pair of bands. A uniform magnetism term Ma^ 
is introduced to shift the bands of H!^ and H'^ upward 
and downward, respectively, without changing the Chern 
number of each band. A gap is opened between two mid- 
dle bands. At half filling, we assume that the Fermi level 
is located in this gap. According to the Chern numbers 
and relative positions of the energy bands, the system 
can be categorized into five cases. Case I: S+ < 2|A| and 
\S-\ < 2|A|. The total Chern numbers of three lower 
bands is 0. However, the nonzero difference between the 
Chern numbers of H!^ and 7J| indicates QSHE. We may 
have two counter-propagating edge states with different 
spins on each edge, although these two edge states do not 
form a time-reversal pair as time-reversal symmetry has 
already been broken. Case II: S+ > 2|A| and S- > 2|A|. 
The Chern numbers of both H!^ and i7| change to zero 
due to the Vs term. When the middle band of is higher 
than the middle band of H'^, the system is in an insulat- 
ing phase as shown in the figure. When the middle band 
of iJ| is lower, the system exhibits QSHE. Case HI and 
Case IV: S+ > 2|A| and l^^l < 2|A|. Case III, when the 
middle band of H!^ is higher than the middle band of H'^, 
the total Chern number is 1. This is a QAHE phase, in 
which there exists one gapless spin-up chiral edge state. 
Case IV: when the middle band of H!^ is lower than the 
middle band of the total Chern number is 0, but it 
gives rise to QSHE. To distinguish Case HI and Case IV, 
one can check the eigenvalues at (7r/2, 7r/2), and we have 
Case HI when Vc + M > 2|A|. Case V: S+ > 2|A| and 
S- < — 2|A|. Due to Vc term, the Chern numbers of both 
H!^ and H'^ change, and the total Chern number is 1. 
Once again the system exhibits QAHE as there is only 
one chiral edge state. 




k k 

X X 

FIG. 3: (Color online) Band structure of Hi; in a strip geom- 
etry. Edge state spectra are shown by dark gray lines (or red 
online), (a) < S+ < 2A. (b) S+ > 2A. 
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FIG. 4: (Color online) Band structure in a strip geometry for 
five different cases. Gray large crosses and dark gray small 
crosses (red and blue online) represent the dispersions for the 
edge states. The columns from left to right are for 
and H!f + H'^, respectively. The horizontal axis is kx G [0,7r] 
and the vertical axis is E/t. 



IV. FERRIMAGNETSIM AND QAHE 

One of the prominent features of the model in Eq.([T]) 
is the appearance of the flat bands due to the unequal 
numbers of sites of the two sublattices, even in the pres- 
ence of spin-orbit coupling. These flat bands give rise to 
the famous flat-band ferromagnetism when the Coulomb 
interaction is turned onM^^ When the system is at half 
filling, the two lower bands are fully filled, and the total 
spin of these two bands is zero since electrons in these two 



bands have opposite spins. The two middle bands are de- 
generate. In this case, if only one single middle band is 
fully filled, the expectation value of the Coulomb interac- 
tion is minimized since the fully polarized electron spin 
in the middle band excludes the double occupancy com- 
pletely at each site. In this way, the ground state of the 
system is ferromagnetic. The total spin is given by the 
degeneracy of the flat band, Stot = Further- 
more, since the antiferromagnetic correlation is dominant 
in the half-filled Hubbard model, this ground state is ac- 
tually ferrimagnetic in which ferromagnetic and antifer- 
romagnetic long range orders coexist*2E When the spin- 
orbit coupling is present, the flat bands will be distorted 
by the ferrimagnetism. When the coupling is strong, the 
ferrimagnetism would significantly distort the fiat bands 
and is suppressed. Let us focus on the case of weak spin- 
orbit coupling, in which the band is expected to be almost 
fiat. It is still possible that the ferrimagnetism could sur- 
vive if the Coulomb interaction is strong enough over the 
band distortion)^ Thus the combination of the flat band 
and the Coulomb interaction provides a reliable mecha- 
nism to realize the spin-dependent staggered potential in 
this system. 

The Hamiltonian with the on-site and nearest-neighbor 
repulsive interactions has the form 

'H^Ho+ Hso + U'^ ni-fn^i+ 

(7) 

^ (n, - l)(n, - 1) - ( 5] - 3^ ), 

where U > is the on-site Coulomb potential, F > is 
the nearest-neighbor repulsive potential, riia = Cia<^ia is 
the number operator for electron with spin a on site i, 
and Hi — Hi-f + Tii^. The chemical potential fj, determines 
the number of electrons in the system. N here is the 
number of unit cells, and the number of electrons is 3iV 
at half filling. In the mean field approximation, the on- 
site interaction is decoupled as'^^ 

where rrii = (clazCi) is the magnetization on site i. Due 
to the asymmetry of the two sublattices, the magnetiza- 
tion on site A and B are different, saying, ttia for the 
sublattice A and tub for the sublattice B. In this way 
the Hubbard term is reduced to 

U nj^n^i = -t/m^y^ ^cr^c — t/mg cJct^q, (9) 

i iGA iGB 

At half filling, since the antiferromagnetic correlation is 
dominant, ruA and tub have opposite signs, based on 
the rigorous results for the Hubbard model 3^1^ This can 
also be illustrated from the calculation of the mean field 
theory. 

The nearest-neighbor interaction may induce the insta- 
bility of charge-density wave (CDW). In the mean field 
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FIG. 5: (Color online) U ~ V phase diagram for A — O.lt 
at half filling. The shadow marks the gapless regions at half 
filling, which should be metallic with strong anomalous Hall 
effect and spin Hall effect. 



approximation, we have 

(10) 

where pA = (rii — l) for the sublattice A and ps = {ni~l) 
for the sublattice B. The CDW order parameter p is given 
by PA = 2p and pb = —p due to the charge conservation. 

After some tedious algebra we may have the zero- 
temperature mean field free energy at half filling 



leads to a set of the mean field equation, 

y ^^i^e(p - i?,(k)) + NUruA = 0, 

i,k 

y ^^^e{fi - E,{k)) + 2NUmB = 0, 
^ omB 

J2 ^^e(M - i5,(k)) + NV{16p - 4) = 0, 

?',k 

J2e{fi- E,(k))-3N = 0. 

i,k 

We solve this set of equations numerically. The calcu- 
lated mean field results are consistent with the rigorous 
results for the Hubbard model. niA and have differ- 
ent signs, which demonstrates the existence of antifcrro- 
magnetic correlation. mA + '2m b ^ demonstrates the 
ferromagnetic correlation. If Wc — 0, tua + 2m b = 1, 
which is one of the rigorous results for the Hubbard 
model. 

Figure. [S] shows the mean field U ~ V phase diagram 
for A = O.lt and t > 0. The CDW order is zero when V 
is small. Vs and M would increase with U and have the 
same sign. The system transits from Case I to Case II 
through a band inversion. However we may see that the 
gap at half filling is not opened in the dashed area of Fig. 
[5l When this gap opens. Case I gives an AFM quantum 
spin hall effect, and Case II is an AFM insulating phase. 
When V is large, the CDW order may become nonzero 
and increase with V. Case HI can be found near the 
transition point where Vc is small, and Case V can be 
found at larger V where Vc becomes large. They both 
present QAHE when the gap opens at half filling. If 
large A is chosen, the Vc term would distort fiat bands so 
much that the orders of Vs and M are suppressed. Thus 
to have a QAHE, we need strong on-site interaction U 
and a nonzero spin-independent staggered field. 



V. CONCLUSIONS 



F{mA,mB, p, fi) 

= Y,^[^'^E,{k)][-^, + E,{k)] + N{Eo + S|l), (11) 

where 9{x) is the step function, £'i(k) is the eigenvalues 
of H{k) in Eq. ^ with H[ = + (v^ + Vc)A + M and 
i7| = Hi — {vs — Uc)A — M, where Vg — U{mA — to_b)/2, 
= ApV, and M = ~U{mA + mB)/2. Eq = Um\/2 + 
Um^ + 8Vp^ — AVp. The summation runs over the whole 
Brillouin zone. The order parameters mA, ms, p, and 
p can be determined self-consistently by minimizing the 
free energy. The variational principle 

6F __ SF _ 5F _ SF _ ^ ^^^^ 
5mA SmB Sp Sp 



We have found that the flat-band ferromagnet may ex- 
hibit QAHE after the inclusion of spin-orbit coupling on 
a 2D decorated lattice. The spin-orbit coupling can in- 
duce topologically non-trivial phase on this lattice, which 
exhibits QSHE. In the present three-band system, the ex- 
istence of the topologically trivial flat band between the 
two non-trivial bands does not affect the formation of 
QSHE. The Coulomb interaction may remove the degen- 
eracy of electrons in the fiat band and lead to sponta- 
neous symmetry breaking, which gives rise to ferrimag- 
netism. The coexistence of ferrimagnetism and CDW 
may break the balance between the helical edge states 
with spin up and spin down in QSHE, and make it pos- 
sible that one branch of edge states is suppressed com- 
pletely, and the other survives. As a result, it gives rise 
to QAHE. 



6 



Acknowledgments HKUST3/CRF/09. 

This work was supported by the Research Grant Coun- 
cil of Hong Kong under Grant Nos. HKU 7037/08P and 



^ N. Nagaosa , J. Sinova, S. Onoda, A. H. MacDonald, and 
N. P. Ong, Rev. Mod. Phys. 82, 1539 (2010). 

2 F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 

^ X. L. Qi, T. L. Hughes, and S. C. Zhang, Phys. Rev. B 78, 
195424 (2008). 

* M. Onoda and N. Nagaosa, Phys. Rev. Lett. 90, 206601 

(2003). 

5 J. Li, R. L. Chu, J. K. Jain, and S. Q. Shen , Phys. Rev. 

Lett. 102, 136806 (2009). 
® C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 

(2005). 

L. Fu and C. L. Kane, Phys. Rev. B 74, 195312 (2006). 

® B. A. Borncvig, T. L. Hughes, and S. C. Zhang, Science 
314 1757 (2006). 

^ Konig, M., S. Wicdmann, C. Briinc, A. Roth, H. Buh- 
mann, L. W. Molenkamp, X. L. Qi, and S. C. Zhang, Sci- 
ence 318, 766(2007). 
^° L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, 
106803 (2007). 

" H. J. Zhang, C. X. Liu, X. L. Qi, X. Dai, Z. Fang, and S. 

C. Zhang, Nature Phys. 5, 438 (2009). 
12 Y. L. Chen, J. G. Analytis, J. H. Chu, Z. K. Liu, S. K. Mo, 

X. L. Qi, H. J. Zhang, D. H. Lu, X. Dai, Z. Fang, S. C. 

Zhang, I. R. Fisher, Z. Hussain, and Z. X. Shen, Science 

325, 178 (2009). 
1^ C. X. Liu, X. L. Qi, X. Dai, Z. Fang, and S. C. Zhang, 

Phys. Rev. Lett. 101, 146802 (2008). 
" R. Yu, W. Zhang, H. J. Zhang, S. C. Zhang, X. Dai, and 

Z. Fang, Science 329, 61 (2010). 

G. Xu, H. Weng, Z. J. Wang, X. Dai, and Z. Fang, Phys. 

Rev. Lett. 107, 186806 (2011). 
1® S. Raghu, X. L. Qi, C. Honcrkamp, and S. C. Zhang, Phys. 

Rev. Lett. 100, 156401 (2008). 
1'' Y. Zhang, Y. Ran, and A. Vishwanath, Phys. Rev. B 79, 



245331 (2009). 

1® K. Sun, H. Yao, E. Fradkin, and S. A. Kivelson, Phys. Rev. 
Lett. 103, 046811 (2009). 

A. Mielke and H. Tasaki, Commun. Math Phys. 158, 341 
(1993). 

2° H. Tasaki, Phy. Rev. Lett. 69, 1608 (1992). 

21 E. H. Lieb, Phys. Rev. Lett. 62, 1201-1204 (1989). 

22 S. Q. Shen, Z. M. Qiu, and G. S. Tian, Phy. Rev. Lett. 72, 

1280 (1994). 

2^ H. Katsura, I. Maruyama, A. Tanaka, and H. Tasaki, EPL, 

91, 57007 (2010). 
2* Z. F. Jiang, R. L. Chu, and S. Q. Shen, Phys. Rev. B 81, 

115322 (2010). 

2^ C. Weeks and M. Franz, Phys. Rev. B 82, 085310 (2010). 
2"^ T. Fukui, Y. Hatsugai and H. Suzuki, J. Phys. Soc. Jpn. 

74, 1674 (2005). 
2^" T. Fukui and Y. Hatsugai, J. Phys. Soc. Jpn. 76, 053702 

(2007). 

2* A. M. Essin and J. E. Moore, Phys. Rev. B 76, 165307 
(2007). 

29 H. Z. Lu, W. Y.Shan, W. Yao, Q. Niu, and S. Q. Shen, 

Phys. Rev. B 81, 115407 (2010). 
^° S. Q. Shen, Inter. J. Mod. Phys. B 12, 709 (1998). 
^1 H. Tasaki, J. Stat. Phys. 84, 535 (1996). 
^2 In the mean field approximation. 



w U {{rii^) (rii^) + Sriii (rii^) + Srii^ (ni^)) 

and the term USuifSnn is droped as a higher order term. 
Define m = (riif) — {riii}. Then (nit) = |(1 + 
(fiii) = 1(1 — "^) when the system is half filled. 



